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Abstract —Most of the available closed-loop explicit guidance algorithms either use constant thrust or 
constant acceleration thrust models. These assumptions which simplify the computation of thrust integrals 
are too restrictive for non-uniform thrusts. The inaccuracy in the thrust modelling works equivalent to 
corresponding perturbations in thrust values, and manifests in large changes in steering angles as well as 
inaccuracy in the injection condition and non-optimality. In this paper, a non-uniform thrust is modelled 
as linearly varying either increasing, decreasing or of constant magnitude over limited time segments. The 
corresponding thrust integrals are developed in a simple way. An earlier developed algorithm is extended 
where each of these linearly varying thrust segment is considered equivalent to a stage of a multi-stage 
rocket while solving the guidance problem. These thrust integrals can be used by most of the explicit 
guidance algorithms and help in improving performance considerably with non-uniform thrust. The 
numerical simulation results show that the algorithms using these thrust integrals provide significantly 
better results as compared to ones using constant thrust models or a model using remaining thrust 
constant. 


t. INTRODUCTION 


In the existing explicit guidance algorithms[l-4] for 
launch vehicles, thrust is assumed to be constant or 
throttled to produce constant acceleration. Often 
these assumptions are too restrictive. In many cases, 
the thrust varies significantly with time, particularly 
in the case of solid fuel rockets. In order to solve for 
the optimum steering problem of closed-loop guid¬ 
ance involving varying thrust, the scheme for con¬ 
stant thrust case can still be used if thrust integrals for 
varying thrust of the form (F, thrust; m, mass; T lo , 
time-to-go; t, time) 


and 



/m(t)} d t 



{F(t)t n /m(t)} df df 


are available. The estimates of these thrust integrals 
are often made based on the average for the stage and 
corresponding mass flow-rate. In a recent study, 
Shrivastava et al.[ 5] used the remaining mass up to 
the end of stage to compute the moving average mass 
flow-rate for the computation of thrust integrals. The 
results corresponding to these two approaches are 
unsatisfactory in accuracy and stability. A varying 
thrust, unless modelled accurately, produces effects 
equivalent to perturbations in thrust and mass and 


(Paper IAF-88-31I presented at the 39th Congress of the 
International Astronautical Federation, Bangalore, India, 
8-15 October 1988. 


hence puts a greater demand on guidance. Conse¬ 
quently this equivalent perturbation, which is felt 
more when the remaining time-to-go is small, results 
in higher injection inaccuracies and instabilities (seen 
as higher turning rate of thrust vector). When only an 
intermediate stage has a varying thrust, but not the 
final stage, the effect is not so critical as compared to 
the case where the varying thrust is in the final stage. 
This necessitates a more accurate estimate of effects 
of thrust. 


2. ANALYSIS 

Any predefined varying thrust can be idealized as 
linearly varying for a limited part of thrust-time 
profile. This way, the whole profile can be divided 
into a number of segments over which the thrust can 
be assumed to vary linearly, (Fig. 1). The number of 
segments depends upon the nature of thrust variation 
and accuracy desired. 

Let the vacuum thrust in the y'th segment, Fj be 
given by 

F{t)=V ex m J (t) (1) 

where V a is constant effective exhaust velocity for the 
entire stage comprising a number of segments, so that 
thrust during y'th segment, F f is directly proportional 
to the mass flow-rate, rhj. Let the absolute value of 
the mass flow-rate for the y'th segment be 

mj = cijt+bj 0<f^7}, (2) 

where a, is the gradient of the mass flow-rate (or 
thrust) during the jth segment, bj is the mass flow-rate 
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Fig. 1. Linear thrust model with four segments for a varying 
thrust stage. 


at r = 0 during the y'th segment, and 7} is the remain¬ 
ing bum time for y'th segment, r is time within y'th 
segment measured from the start of the segment (if 
the ignition of segment j has not started) or from the 
current time (if the current segment is y'th segment). 
The mass at any time r during the y'th segment is given 
by 

m j(t) = Cj - 1 12ajt 2 - bjt (3) 

where c y is the mass of the vehicle at t = 0 during the 
y'th segment. Since c y and m y represent mass, 


nij > 0 . 

The thrust in the segment j at anytime t is obtained 
by substituting eqn (2) into (1) and is given by 

Fj{t) = V a (a jt + bj). (4) 

Evaluations of a Jt b jt Cj are done a priori by making 
a linear least squares fit of the thrust over each 
segment. Since V a is known for a given motor, a, and 
bj can be computed from 

Fj(t) = Tjt + Sj 

where r y r + i y is the line representing thrust in the 
least squares way in y'th segment. Then, using eqn (4) 
one gets 


and 

a, = r jl V ao 

(5a) 


II 

X 

(5b) 

Evaluation of thrust integrals 


The thrust integrals for the y'th segment 
as follows: 

are defined 


\ Tj {Fj(t)rimj(t)} dr 

to 

(6a) 

^Drtj 

Tj f ! {Fj(t)C/mj(t)} dr dr, 

0 Jo 

(6b) 

where 7} is the burn time for segment y. 



As the mass flow-rate has been considered to be 
positive, bj > 0. On the other hand, a, may be posi¬ 
tive, negative or zero in a segment depending on the 
nature of variation of thrust. When a y = 0, it repre¬ 
sents a constant thrust over the segment and the 
relations obtained for thrust integrals for constant 
thurst[l] are valid. For the cases of a, > 0 and a y < 0, 
an analytical solution for single integrals, I s „., are 
derived next. The double integrals, I Dnj , can be ex¬ 
pressed in terms of I Snj by carrying out integration, by 
parts, of eqn (6b) and using eqn (6a) as 


= T,L 


j Snj 




Let 


x —bj/Oj + r (7a) 

and 


l 2 —bj/aj + Icj/Oj. (7b) 

Case 1. Linearly increasing thrust (a y > 0) 

On substituting expressions for Fj, ntj in expression 
for /$,, given by eqn (6a), one gets 

Isn t = V„ r fiajt + bj)l(Cj - \l2ojt 2 - bjt) d t. (8) 

Rewriting: 

t n /(Cj - l/2«/ 2 -V) = -2 r-'/aj + {2(c y /a y )f" -2 

- 2(h y /fl y )r"- *}/(c y - l/2a y f 2 - bjt). 

On substituting the above relation in eqn (8) one 
obtains 


2T 2bjt"~' Jt 

n Oj(n - l)Jo 


+ 2(h y /a y )/j(„_i) y — 2(Cj/a j )I S ( n _2 jj - (9) 

^Snj can be solved recursively for n ^ 2 if and 4i y 
are known. Equation (8) can be rewritten as 


t" d/dr{ln(c y — l/2a ; / 2 — bjt)} dr. (10) 


For a y > 0, 1 2 > 0 as c y > 0. So, / is a real number 
which is taken to be positive square-root of 
(bj/aj + 2Cjlaj). On using eqns (7a) and (7b), one 
obtains 


m y (f) = Cj — 1/2 Ojt 2 — bjt= Oj(l 2 — x 2 )/2 >0 (11) 

which leads to (/ — x) > 0. 

From eqn (7a) one gets 

r" = x"[l - bjKajx)]" 

1+1 02 ) 

where 



n Cr = n!/{(n — r )!r!}. 
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On substituting eqns (11) and (12) into eqn (10) 
and expanding the terms one gets 

-4. > /^ex = | x " d/d/ {ln(/ 2 — x 2 )} dt 
— nbj/aj J x"~ 1 d/dt {ln(/ 2 — x 2 )} dt 


+ ■■• + ( 
Let 


-i rwJJ ' 


d/d/{ln(/ 2 — jc 2 )} dr. (13) 


— J*^ x " d/dt {ln(/ 2 — x 2 )} dt. 

Equation (13) can be rewritten as 

hrtj = — ^exI/lSny ~ n bj! a jI'lS(n - (Jj + ’ ' ' 

+(-mb J /a J yr lS0j ]. (14) 

Equation (14) can be evaluated for ! S0j and Isij by 
substituting n =0 and 1, respectively, as follows 

= — VcxI'lSOj 

hi/ ~ ~ ^exOlSly — bjl’isu)j/(lj). 

After integrating by parts for integral I\ Snj one 
obtains 

— « 1 ln(/ 2 — jc 2 )d/. (15) 

On using eqn (15), one can evaluate I' lSn/ for n = 0 
and 1 as follows 


hsoj = fln(/ 2 - x l )]l‘ 

hsij = [* ln(/ 2 - x 2 )# - [x ln(/ 2 - x 2 ) - 2x 

4- / ln{(/ + x)/(/ - x)}]f 
where x and / are defined by eqn (7). 

Case 2. Linearly decreasing thrust (a, < 0) and l 2 > 0 

Since / 2 >0, / is a real number and is taken as 
positive square root of l 2 . w,(f) can be written as 

mj(t) = (— l/2fly)(x 2 — l 2 ), (16) 

for Oj < 0, (x 2 — l 2 ) > 0. 

So, (x — /) and (x + 1) are both of the same sign, 
either positive or negative, and hence 
ln{(x — l)/(x + /)} is always defined. 

On using eqn (16), a relation similar to eqn (13) can 
be written as 

+ i »c,(- l) r {tya,*} r ] d/df(x 2 - 1 2 ) d t. (17) 
Let 


where integrals r 2S n, correspond to Case 2. Equation 
(17) can be rewritten as 


— — VtxU'lSn/ — . l)j 

+•••+(-(is) 

for n ^ 1 is obtained from eqn (17) as follows: 

hi, = — Vtxl’isoj 

htj = ~ ^exKiSlj — bjhsVjl 0 /]!'■ 

On integrating the expression for /J^, by parts, one 
gets 

I’lsnj = [*" ln(x 2 - - n j* ln|(/ 2 - x 2 )|x" “ 1 At. 

(19) 


Using eqn (19) one finds 


I'isnj = [x ln(x 2 - 1 2 ) ]j> 

hsij = [x ln(x 2 - / 2 )]f - [x ln|(/ 2 - x 2 )| 

Tj 


where x and / are defined by eqn (7). 


-2 x+n 


l+x 


Case 3. a y <0 and l 2 <0 
Let 

l^^-ibj/aj + lcj/a/) (20) 

then 

l' 2 > 0 for (b 2 + 2 C/Oj) < 0 
so 

irij{t) = cj - 1/2 Ojt 2 - bjt = (-1/2o,)(x 2 + 1' 2 ). 

On using eqn (20), a relation similar to eqn (13) can 
be written as 



+ Z «q(~ mbj/ajxY} d/dt {ln(x 2 - 1' 2 )} dt. (21) 

r« 1 

Let 


hsn t — x " d/dt {ln(x 2 4- 1' 2 )} dt 

where corresponds to Case 3. 

Then 


Isa/ ~ ~ VtxV'lSn, — n bj/ a jI')S(n-l)j + ' ' ' 

+ (-mbj/ajyr 3SOj }. (22) 

V can be solved for n >2 if I x , I slj are known. After 
substituting n = 0 and 1 in eqn (22), one gets /*, for 
n < 1 as follows 


/« 


j> 


d/dt{ln(x 2 -/ 2 )} dt 


IsOj — ~ ^cxhsO/ 

h)/ = — ^ext^JSIy —b)! a ll 3S0y]- 
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Integration by parts of the expression for I' 3Shj yields, 

1’lSn, = [*" ln<y + 1 2 )$ 

~ j Tj nx "~ 1 ln(x 2 + 1' 2 ) d/. (23) 

The integrals given by eqn (23) can be evaluated for 
n = 0 and 1 as follows: 

-/*>]? 

I'isij = [* M* 2 + l' 2 )]o J - lx In(jc 2 + 1 2 ) 

— 2x + 21' tan _l (x//')]o J where x is defined by eqn 
(7) and /' is defined by eqn (20). 

3. SOLUTION OF THE GUIDANCE PROBLEM USING 
VARIABLE THRUST 

Having obtained the thrust integrals for each seg¬ 
ment the system can be treated as a multi-stage rocket 
where each segment works like an equivalent single 
stage. The solution where guidance is operative in 
more than one stage is treated in detail in [1]. Using 
the method developed, changes in position and veloc¬ 
ity can be determined with the number of segments 
in the varying thrust stage taken to be equivalent to 
number of stages. 

Using the algorithm for the solution of guidance 
problem given in [1], effects of thrust in terms of 
changes in velocity and position corresponding to 
each of the segment for a given set of guidance 
parameters can be determined. The corresponding 
changes due to gravity can also be estimated using 
Encke’s method[l]. 

Effects of thrust 

The required change in velocity to be inparted by 
thrust is given by 

VthR 4 V r -V 0 -V g 

where v T is the required velocity at injection, v 0 is the 
current velocity and 



v g can be estimated using Encke’s method[l]. (Bold 
letters indicate vectors.) 

The optimal unit thrust vector, developed in {1] is 
given in terms of unit vectors I K and l yc by 

■(0 “If*+ !*«(* — t c )]/[1 +pj(t- i c ) 2 ] l/2 (24) 

where and I*, form the plane of correction, is the 
tangent of the steering angle made from I xc , and t c is 
a time constant. The denominator of eqn (24) can be 
written as 

[i+pl(t-T3 t ] ll2 =i-il2ptit-Z'Y 

+ 3/8/>e(r — to) 4 + ''' ■ (25) 
The estimated change in velocity imparted by thrust 


for a specified set of guidance parameters for a 
/^-segment varying thrust stage is given by 

VthA 4 L Z f ^/K{1 +p](t - t C; ) 2 } 1/2 ] dt 

j - c.s. J 0 

where c.s. is current segment, is the direction in 
which velocity constraint is met, and 7) is the remain¬ 
ing bum time for the yth segment. x Cj for the current 
segment can be computed using the following relation 

0 = Z f F jP»0 - t c,)/K{1 +Pe(t - T Cy ) 2 } 1/2 ] df. 

y*c.s. Jo 

x Cj of any other segment can be computed with respect 
to t Cj of the current segment as follows 

T c J+ , = tcj — Atj + „, 

where x c . is for the current segment and A t j+n is the 
time-difference between ignition time of (j + n) seg¬ 
ment and the current time. 

The expression for r lhA can also be written in a 
similar manner on using the following relation 

r . hA=z r r f 

>-c.s. LJo Jo 

x {1 +Pe(t -T c> } 1,2 ]dsd/ 

+ b, J F j(t)/[mj(t){l +pl(t - x Cy ) 2 } ,/2 ] di 

+ K Z [ f [ F j( s )Pe( s ~ t c )/H(s) 

LJo Jo 

x {1 +/^(s -r C; ) 2 } l/2 ]ds At 

+ b, | F ,(t)l[mj(t){ 1 +p](t - t c> ) 2 } 1/2 ] d t 

where t r/ is the remaining time for the thrust cut-off 
after the end of Jth segment. The corresponding 
required change in position to meet the end con¬ 
straints can be given as 


UhR r o Vo F go tg, 

where r T is the position at injection, r 0 is the current 
position, and r g is defined as 

g{r(s)} ds d t. 

r, can be estimated using Encke’s method[l]. 

For the solution, the following constraints must be 
satisfied 



and 


v thR — v thA> 


r thR — r thA* 

Using eqn (25), the expression for and can 
be written in terms of l m . and I Dnj for each segment 
and thus v,^ and can be evaluated. 

The steering parameters used in eqn (24) can be 
determined after forming algebraic constraint equa¬ 
tions derived from vector constraint equation given 
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Table 1. Variable thrust-time and mass 
profile 


Time from 
guidance 
initiation 
(s) 

Mass 

(kg) 

Thrust 

(N) 

0 

4300.00 

11250 

25 

4201.04 

12500 

50 

4092.71 

13500 

75 

3975.63 

14600 

100 

3846.67 

15650 

125 

3712.71 

16500 

150 

3572.09 

17250 

175 

3425.52 

18000 

200 

3282.50 

18600 

225 

3125.83 

19000 

250 

2966.67 

19000 

275 

2812.51 

18000 

300 

2668.76 

16500 

325 

2538.34 

14800 

350 

2420.00 

13400 

365 

2367.67 

9800 

375 

2330.75 

8300 

400 

2300.00 

0 


above. These constraint equations are solved using 
differential correction method[l] where partial deriva¬ 
tives are found analytically. 

4. NUMERICAL SIMULATIONS AND RESULTS 

In order to investigate the effectiveness of the 
algorithm with varying thrust, numerical simulations 
are conducted. For this, the performance of the 
guidance in the final stage with varying thrust and 
having thrust cut-off facility is considered. 

The variation of mass of the vehicle and thrust¬ 
time profile of this stage are given in Table 1. Figure 
2 shows the variation of the thrust with time consid¬ 
ered for simulation. Table 2 gives the initial position 
and velocity of the rocket at the start of guidance. 
The desired orbit at injection is also given in Table 2, 
which corresponds to a 900 km circular Sun-syn¬ 
chronous orbit. 

Table 3(a) shows the idealization of the varying 
thrust given in Table 1 by different models. Three 

Linear varying 



Fig. 2. Thrust-time profiles using a linear time-varying 
thrust model, constant thrust model and moving average 
constant thrust model. 


Table 2. Initial and target conditions of the vehicle 


Initial Position 

(2181127.34,-1978553.22, 

(x,y,z components) 

6570848.69) m 

Initial velocity 

(3519.09498, -4212.32, 

(x,y,z components) 

-1766.37778) m/s 

Geocentric distance at injection 

7278137.0 m 

Inertial velocity at injection 

7440.46107 m/s 

Flight path angle at injection 

0.0° 

Inclination at injection 

98.98” 

Longitude of the ascending node 
(measured from the Greenwich 
meridian at lift off) 

263.76” 


such models are considered. The first model (model 
A) corresponds to the linearly varying thrust model 
developed in this paper, where the thrust-time profile 
is divided into three segments. The first segment is 
from the start of the state up to 250 s where thrust is 
assumed to vary linearly from 12,000 to 20,000 N. In 
the next segment, from 250 to 350 s, the thrust is 
assumed to vary from 20,000 N to 13,500 N linearly. 
In the last segment, from 350 to 400 s, the thrust is 
assumed to drop from 13,500 N to zero linearly. The 
corresponding values of coefficients a,b,c for these 
three segments are shown in Table 3(b). 

The model is investigated for both increasing and 
decreasing thrust. The other two models considered 
for comparing its performance are similar to those 
discussed by Shrivastava et al.[S\. Model B corre¬ 
sponds to a constant thrust of 15,000 N for the entire 
stage which has nearly the same total impulse as the 
varying thrust given in Table 1; and model C corre¬ 
sponds to a constant mass flow-rate obtained by 
taking the ratio of remaining propellant in the stage 
and the maximum possible duration of bum of the 
stage. For this duration the thrust, which is computed 
corresponding to the moving average mass flow-rate, 
is assumed to be constant. Figure 2 also shows the 
thrust-time profile assumed by model B and the 
constant thrust considered by Model C at 200 s after 
the ignition of the stage which has a maximum 
duration of bum equal to 400 s. 

For investigating the performance of these models, 
the guidance algorithm described in [1] is used. The 
first case corresponds to an ideal situation where the 


Tabic 3. (a) Assumed thrust-time profiles by models A,B,C 


Time from 
guidance 
initiation 
(s) 

Thrust 
using 
model A 
(N) 

Thrust 
using 
mode) B 
(N) 

Thrust 
using 
model C 
(N) 

0 

12000 

15000 

15000 

200 

18400 

15000 

14737.5 

250 

20000 

15000 

13333.2 

350 

13500 

15000 

7200 

400 

0 

0 

0 


Table 3. (b) Coefficients a,b,c for three segments for Model A 


Time from 


guidance 

initiation 

(s) 

Segment 

No. 

0) 

Coefficient 

a J, 

(kg/s 2 ) 

Coefficient 

b j 

(kg/s) 

Coefficient 

c ) 

(kg) 

0 

1 

0.01066667 

4.0 

4300 

250 

2 

-0.023333 

6.66667 

1966.667 

350 

3 

-0.093333 

4.66667 

2416.667 






Pitch angte (deg.) 
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Table 4. Injection errors using different models 



Ideal 

Model A 

Model B 

Model C 

Bum time (s) 

366.21 

366.03 

366.23 

366.954 

Error in semi-major axis (m) 

-0.009 

-2.12 

16561.9 

497.6 

Error in eccentricity 

2.62 x 10“ 8 

6.0 x 10 _! 

2.52 x 10~ 3 

1.10 x 10“ 3 

Error in injection position 

(0.075,-0.12,-0.13) 

(177.22, -274.50, 

(-4497.41, 7081.09, 

(3214.63,-5084.77, 

(; K,y,z component) (m) 


291.90) 

7529.4) 

-5422.06) 

Error in injection velocity 
(x,y,z components) (m) 

(-0.11,0.16,-9.1) x 10“ 6 

(-0.20,0.82,1.04) x 10~ 3 

(1.53,-9.69,-2.63) 

(-0.07,0.20,-0.21) 


actual thrust is identical to the linear thrust of model 
B. The injection accuracy and bum time for this case 
are given in Table 4. The table also shows the 
injection accuracies and bum times when thrust 
corresponding to models A, B and C, respectively, are 
used. In the ideal case, the injection error is nearly 
zero. The errors when model A is used are 2.12 m in 
a semi-major axis and 6.0 x 10' 5 in eccentricity which 
are quite small. As compared to these results, when 
models B and C are used, the errors in semi-major 
axis are 16561.9 and 467.6 m, respectively. The errors 
in eccentricity using models B and C are 2.522 x 10 3 
and 1.103 x 10 -3 , respectively. The pitch steering 
angles using the ideal case and models A, B and C are 




Fig. 3. Pitch (a) and yaw (b) angle profile using actual thrust 
identical to model thrust, linear time varying constant thrust 
and moving average constant thrust models. 


given in Fig. 3(a). Figure 3(b) shows the correspond¬ 
ing plots for yaw angles. It can be seen that the ideal 
case and model A follow almost linear variation of 
steering angles, whereas in case of B and C, the 
instability sets in after about 270 s of guidance and 
the situation worsens thereafter. To use models B and 
C without instability one has to remove position 
constraint nearly 100 s before the actual cut-off. This 
will introduce large injection errors. Thus, for a case 
where the final stage has a variable thrust, model A 
alone provides a satisfactory solution. 

5. CONCLUSION 

In the published literature, the closed-loop guid¬ 
ance algorithm has been developed for constant 
acceleration and constant thrust for a rocket. In this 
paper an existing algorithm has been extended to 
include a varying thrust case. A model consisting of 
segments over which thrust varies linearly is used. 
The corresponding thrust integrals are developed for 
linearly increasing and decreasing thrust. This can be 
used along with most of the guidance algorithms 
developed earlier. The numerical simulation results 
show that the thrust integrals developed for varying 
thrust give significantly better results as compared to 
those using constant thrust model and the model 
using remaining thrust constant. 
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